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Abstract 

We use an elastic rod model with contact to study the extension versus rotation diagrams 
of single supercoiled DNA molecules. We reproduce quantitatively the supercoiling response 
of overtwisted DNA and, using experimental data, we get an estimation of the effective 
supercoiling radius and of the twist rigidity of B-DNA. We find that unlike the bending 
rigidity, the twist rigidity of DNA seems to vary widely with the nature and concentration 
of the salt buffer in which it is immerged. 
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At first the DNA molecule can be seen as the passive carrier of our genetic code. But in order 
to understand how a 2 m long string of DNA can fit into a 10 fim nucleus, it is clear that one has 
to consider the mechanical properties of the molecule. Namely the fact that the DNA double helix 
is a long and thin elastic filament that can wrap around itself or other structures. Furthermore, 
the elastic properties of DNA play an important role in the structural dynamics of cellular process 
such as replication and transcription. Although the basic features of DNA were elucidated in the 
years following the discovery of the double helix geometry, it is only during the past decade that 
few groups, using different micro-techniques, have been able to manipulate single DNA molecules 
in order to test their mechanical properties. 

A first way to manipulate a single molecule of DNA is to attach a polystyrene bead at each 
of its two ends. One bead is then stuck to a micropipette and the other is held in an optical trap. 
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One pulls on the molecule by moving the pipette and measures the force through the displacement 
in the trap. This artefact is called an optical tweezer pQ. Of course elastic properties of a DNA 
molecule will in general depend on the sequence of base pairs (bp) it is made of. Nevertheless for 
long molecules, i.e. more than a hundred bp, the behavior of the molecule can be described by 
a coarse-grained model known as the worm like chain In this model, DNA is considered as a 
semi-flexible polymer with bending persistance length A. This is the contour length over which 
correlations between the orientation of two polymer segments is lost. It can be viewed as the ratio 
of the elastic bending rigidity K to the thermal energy fcgT. The common accepted value is 
A = 50 nm in physiological buffer. 

Another way of manipulating a single DNA molecule is to use a magnetic tweezer [Sj- Here 
the molecule is locked on a glass surface at one end, and glued to a magnetic bead at the other 
end. The bead is now controlled by a magnet that one rotates in order to input a twist constraint 
in the system. The pulling force is tuned via the monitored distance between the magnet and 
the bead. The end-to-end distance of the DNA molecule is measured thanks to a microscope. 
Experiments are carried under constant pulling force /. One gradually rotates the magnet around 
an axis perpendicular to the glass surface while monitoring the relative extension of the molecule 
z = Z/L where L is the total contour length of the molecule and Z is the distance in between the 
bead and the glass surface, see fig. (0). Force is measured using the Brownian motion of the bead. 
No direct measurement of the twist moment is possible with magnetic tweezers. Consequently the 
twist persistence length C is not directly available. 




n = n>0 

Figure 1: The magnetic tweezer experiment. 

When no rotation is put in, the DNA molecule behaves like a semi-flexible polymer, i.e. the 
relative extension 2 is a function of the temperature T, the persistence length A and the applied 



2 



pulling force /: 

'Ak B T\^ 



*(n = 0) = l-(^-J . (1) 

This relation is used to extract the bending persistence length A from experimental data. Then 
under gradually increased rotation, the extension z decreases with the number of turns n put 
in and eventually the molecule starts to wrap around itself. Geometrically speaking, the DNA 
molecule is coiling around itself in a helical way. Since the molecule is already a double helix, we 
speak about supercoiling. Each helical wave of the super helix is called a plectoneme. Different 
attempts has been made to model this experiment, introducing the concepts of wormlike rod 
chain £TJ. or torsional directed random walk [5] but without taking into account the possibility of 
contact, or dealing with the plectonemes in an ad hoc way j^j. 

Here we present an elastic model that specifically includes the self-contact of the DNA 
molecule, but does not take into account thermal fluctuations. Our point is that, in the regime 
where plectonemes are formed, the relevant physical information is already present in our zero- 
temperature elastic rod model with hard- wall contact. Due to the base pair conformation, it seems 
natural to consider DNA as an elastic rod with a non-symmetric cross-section (i.e. two different 
bending rigidities). Nevertheless it has been shown that when dealing with a long enough molecule 
(more than a few dozen base pairs), one could simplify the model and work with an effective rod 
having a symmetric cross-section (i.e. one bending rigidity that is the harmonic mean of the two 
raw rigidities) [7j. In order to focus on supercoiling, we consider the simplest elastic rod model 
than includes twist effects and can have 3D shapes. Following the classic terminology of Euler 
planar elastica for twistless 2D shapes, we call the present model the Kirhhoff ideal elastica [H|. 
The elastic energy reads: 



E= \ I (K K 2 (s) + K 3 T 2 (s))ds 
z Jo 



where s is the arc-length, k(s) the curvature of the centre line, r(s) the twist rate of the cross 
section around the centre line (which in the case of an ideal elastica is a constant of s), and K and 
K 3 are the bending and twist rigidities respectively. The Kirchhoff equilibrium equations read: 

F'(s)+p(s) = (2) 
M'(s) + R'(s) x F(s) = (3) 

where F(s) and M(s) are the internal force and moment respectively. The external force per unit 
length p(s) can model an electrostatic repulsion, gravity or hard-wall contact. The centre line of 
the rod is given by R(s) and t(s) = R'(s) is its tangent. In the case of an ideal rod it can be 
shown |Hj that: 

K t'(s) = M(s)xt(s) (4) 
K Q dy'(s) = (M(s)-r(K 3 -K )t(s))xd 1 (s) (5) 

where di(s) is a unit vector, lying in the cross section, that follows the twist of the centre line. 
For a DNA molecule, it is generally taken as the vector pointing toward the major groove. For 
the parts free of contact, we have p(s) = 0. In case of self-contact there are two points along the 
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z = Z/L 




Figure 2: Fitting the linear part of the response curve. 



rod, say at arclength si and s 2 , where the inter-strand distance |-R(si) — R(s2)\ is equal to twice 
the radius of the circular cross-section (which we denote by p). At point s\, we introduce a finite 
jump in the force vector F(s): 

F{s < Sl ) = F(s > Sl ) + ^ (R( Sl ) - R(s 2 )) (6) 

where AF 12 is a positive real number. The same treatment is done at point s 2 , with the same AF 12 
[HUH!]- This corresponds to having a Dirac function for p(s) in (J2J). In case of continuous sections 
of contact p(s) is a function with changing direction and intensity. In our model we only consider 
cases where the contact occurs either at discrete points or along straight lines. This model has 
already been used [H1H2] an d is weu described in [TUj . 

Ideally we are looking for equilibrium configurations of the ideal elastica that match the 
boundary conditions corresponding to the magnetic tweezer experiment. Nevertheless, in case of 
a long rod (L ^> p) with a large number of turns n put in, the exact geometry of the boundary 
conditions are less important, and we choose clamped ends as used before We numerically 
find equilibrium configurations matching the boundary conditions using classical path following 
techniques; first a self made algorithm relying on multiple (or parallel) shooting, then using the 
code AUTO [T3] that discretises the boundary value problem with Gauss collocation points. The 
different types of solutions (straight, buckled, supercoiled) are found in the following way. We fix 
the radius p and the vertical component / of the force vector F(s = 0) acting on the bead. We 
first consider a straight rod with no rotation n = 0. Then we twist the rod by gradually increasing 
n. At a certain bifurcation value, the path of straight solutions crosses another path of buckled 
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solutions. Following this new path, we end up crossing another path of solutions with one discrete 
contact point. This latter path will itself intersect a path of configurations with two contact 
points. Solutions with up to three discrete contact points have been found p]. Eventually the 
configurations with three discrete contact points bifurcate to solutions including a line of contact 
in addition of discrete contact points. We call them supercoiled configurations. In a supercoiled 
configuration, the parts of the rod which are in continuous contact have an helical shape. We 
call this twin super-helix a ply. The super-helix is defined by its radius p and its helical angle 9 
(see fig. (0)). Each time we change the fixed pulling force / or the radius p, the entire numerical 
continuation has to be rerun. Since we do not consider thermal fluctuations, the first part of our 
numerical response curve does not correspond to what is found experimentally. On the other hand 
our model reproduces quite precisely the part of the response curve where the distance z decreases 
linearly with the number of turns n, provided we identify p not with the crystallographic radius of 
the DNA molecule but with an effective supercoiling radius due to electrostatic as well as entropic 
repulsion. We numerically find that in the linear regime the helical angle 9 does not vary with n. 
We fit numerical solutions curves as in [14J and we find that 9 only depends on /, K , and p: 

f = ^ fo{9) with <f) 3 (9) = 1.65805 9* . (7) 

This result enables us to extract the effective supercoiling radius p and the twist rigidity K 3 
from magnetic tweezer experiments on DNA. First we note that the number of turns n applied 
to the magnetic bead can be interpreted as the excess link of the DNA molecule: n = ALk. Link 
is normally defined for a closed ribbon but carefull use of a closure permits the introduction of 
the link of an open DNA molecule [751 ITS'] . We use the Calugareanu- White- Fuller theorem to 
decompose the excess link: 

(n =) ALk = ATw + Wr (8) 

where ALk (resp. ATw) is the difference between the actual link (resp. twist) and the intrinsic 
link (resp. twist) of the double helix. The writhe Wr is the average number of crossings of the 
centre line one sees when looking at the molecule from (all the possible) different viewpoints. In 
the ply part, the twist rate is related to the helical angle 9 via a mechanical balance [TTH ITTj: 

A " (tan 29 -sin 29) , (9) 



2pK, 



where e = ±1 stands for the chirality of the ply ^7]. Since the twist rate is constant along the 
rod, we have ATw = tL/(2tt). Generally writhe is not additive but using Fuller theorem ^H] 
with a carefully choosen reference curve we may write : 

Wr = Wr loop + Wr taUs + Wr ply . (10) 

Here we neglect Wri oop and Wr ta u s . Carefully choosing a reference curve to apply Fuller theorem 
or directly computing the writhe from the double integral yields |15j : 

Wr ply = (-e)^-sm29. (11) 

47T p 
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Figure 3: Experimental response curves with a 48kbp DNA molecule in a 10 mM phosphate buffer 
Each curve corresponds to an experiment carried at a constant pulling force. Forces are, from 
bottom to top, 0.25, 0.33, 0.44, 0.57, 0.74, 1.1, 1.31, 2.2, 2.95, 4.15, and 4.9 pN. 
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The total contour length (or number of base pairs) L of the DNA molecule is given and we write: 



Lpiy L Li 00 p Lt a n s 
We neglect Li oop and we use an heuristic equation for L taas : 

L 



-'tails 



Z a) 



z(oy 



(12) 



(13) 



which means that when the molecule is supercoiled the tails parts, which are not supercoiled, are 
as disordered as the molecule was when no rotation was put in. For positive supercoiling [n > 0) 
we get a left handed ply (e = —1). From (jSJ), ©, and (fTTj) we obtain: 



ALk 



^ (tan 26 



sin 26) + sin 26 1 



Z(0) 



(14) 



Introducing the supercoiling ratio a = ALk/Lko and taking h = 10.5 base pairs per turn and 
Sz = 0.34 nm of rise for each base pair, we get an helical pitch for the DNA double helix H — h8z = 
3.57 nm and we have: Lko = L/H. Equation (|T4^) can then be inverted to give an approximation 
of the linear part of the response curve in the (a, z) plane : 



z(0) 



1 + 



Kz V cos 26 



1 



4np 



Hsin26 



(15) 



We use (J7J) and (|15|) to extract information from experimental response curves of extension- 



/(pN) 


6 (rad) 


p(nm) 


K 3 /K 


0.25 


0.427 


6.85 


1.88 
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0.449 


6.60 


1.87 
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0.467 


6.17 


1.92 
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5.47 


1.84 


0.74 


0.504 


5.55 


2.01 


1.1 


0.488 


4.26 


1.86 


1.31 


0.471 


3.64 


1.58 


2.2 


0.503 


3.21 


1.65 


2.95 


0.507 


2.81 


1.62 



Table 1: Results for the 48kbp DNA molecule in a 10 mM phosphate buffer. 

rotation experiments. Since experiments are carried at given (fixed) / and since K is obtained 
from Q with K Q = AksT, then from each experimental curve in the (a, z) plane we only need 
: (i) the relative extension at o = 0, which we denote by z(0), (ii) the slope of the linear part, 
which we denote by a, (Hi) the ordinate at the origin of the straight line fitting the linear part of 
the response curve, which we denote by zp. Using and (JTHJ) we get an equation for 6 : 

^z(0) Ik 

a = -H^2e\T Me) (16) 
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Figure 4: Experimental response curves with a llkbp DNA molecule in a 150 mM phosphate 
and 5 mM magnesium (Mg2+) buffer. Forces are, from bottom to top, 0.45, 1.45, and 4.3 pN. 
Unpublished data from Gilles Charvin. 



8 



Once we know 6, we can extract the effective supercoiling radius p from (JJJ) and the effective 
stifnesses ratio from ([T5]): 

*> = (-}—l)/(jL-l) (17) 



/(pN) 


9 (rad) 


p (nm) 


K 3 /K 


0.45 


0.307 


2.79 


1.09 


1.45 


0.319 


1.67 


1.00 


4.3 


0.348 


1.15 


0.99 



Table 2: Results for the llkbp DNA molecule in a 150 mM phosphate and 5 mM magnesium 
(Mg2+) buffer. 

Results are shown in table ^ and 121 and can be compared to results of JH] where another 
micro-technique was used and a value of K 3 ~ 100 nm k%T (in a lOOmM NaCl + 40 mM Tris-HCl 
buffer) was found. From our results, we see that the effective supercoiling radius decreases with 
the intensity of the pulling force, and with the strength of the salt buffer. Fitting experimental 
data with (0) yields K = 51 nmksT for fig. Eland K = 57 nm ksT for fig. HJ This means 
that the bending rigidity does not vary widely with the salt properties of the buffer. On the 
other hand, present results imply that the twist rigidity strongly decreases with the salt strength. 
This in turns implies that the electrostatic contribution (repulsion of the charges lying on the 
sugar-phosphate backbone) to this effective twist rigidity is rather important. 

It is a pleasure to thank G. Charvin, V. Croquette, and D. Bensimon for supplying us with 
data from experiments. 
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